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^ ; ABSTRACT 

First using a sample of 32 GRBs with known redshift (Guidorzi et al. 2005) and then a sam- 
ple of 551 BATSE GRBs with derived pseudo-redshift (Guidorzi 2005), the time variabil- 
' ity/peak luminosity correlation {V vs. L), originally found by Reichart et al. (2001) using a 

^ : sample of 18 GRBs, was tested. For both samples the correlation is still found but less rel- 

evant due to a much higher spread of the data. Assuming a straight line in the log L-log V 
plane (logi — mlogV^ + h), as done by Reichart et al., both Guidorzi et al. and Guidorzi 
found that the line slope for both samples is much lower than that derived by Reichart et 
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^ ■ al.: m I-HqI (Guidorzi et al. 2005), m = 0.85 ± 0.02 (Guidorzi 2005), m = S.Stg g 

, (Reichart et al. 2001). Reichart & Nysewander (2005) discuss our results and attribute the dif- 

ferent slope to the fact we do not take into account in the fit the variance of the sample (also 
called slop), and demonstrate that, using the method presented by Reichart (2001), the ex- 
panded data set of Guidorzi et al. (2005) in log L-log V plane is still well described by a line 
Q ■ with slope m — 3.4lg g. Here we compare the results of two methods accounting for the slop 

|H \ of the sample, the method implemented by Reichart (2001) and that by D' Agostini (2005). We 

C/3 . demonstrate that the method used by Reichart et al. (2001) to estimate the straight line slope, 

' provides an inconsistent estimate of the parameter when the sample variance is comparable 

with the interval of values covered by the GRB variability. We also demonstrate that, using 
. . the D' Agostini method, the slope of the logL-log correlation is still consistent with that 

' derived by us earlier and inconsistent with that derived by Reichart & Nysewander (2005). 

^ , Finally we discuss the implications on the interpretations proposed for the V — L correlation 

and show that our results are in agreement with the peak energy/variability correlation found 
by Lloyd-Ronning & Ramirez-Ruiz (2002) and the peak energy/peak luminosity correlation 
(Yonetoku et al. 2004; Ghirlanda et al. 2005). 
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1 INTRODUCTION 

Most of our knowledge about the Gamma-Ray Burst (GRB) phe- 
nomenon is derived from their spectra and light curve profiles, but 
it is recognised that other observational probes (e.g., polarisation of 
the gamma-rays) would give key information to the solution of the 
GRB enigma. Among these probes, it is recognised the importance 
of the erratic time variability of the GRB time profiles. For exam- 
ple, in the GRB internal shock model, not very variable radiation is 
expected to be produced at radii lower than the p hotosphere radius 
in wh ich the shocks remain optically thin to pairs iKobavashi et alj 
l2002h . while highly variable radia tion is expecte d to be produced 
in shocks above this radius ilMeszaros et all2002h . Also in the sub- 
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jet model by loka & Nakamura time variability is expected 

and its amplitude related to the viewing angle of the burst. 

A key objective of the study of GRBs is to establish whether 
GRBs can be reliably used as standard candles and to determine 
the optimal relationship between observed and intrinsic properties. 
The recent discovery of GRBs with bright afterglows at redshifts 
z > Q highlights their power as probes of the high redshift Universe 
(Haislio et al. 2006; Kawai et al. 2006), but spectroscopically- 
confirmed redshifts are only a fraction. In contrast, the characteris- 
tics of gamma-ray light curves, which are available for all GRBs, 
offer a potentially independent estimate of luminosity distance for 
statistically-significant samples. 

An important problem, however, in deriving the intrinsic cor- 
relations is the interpretation of scatter in correlations, which may 
be produced by measurement methods, construction of samples 
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with properties measured by satellites with differing response func- 
tions, small samples and different statistical analysis methods and 
intrinsic physical differences in the GRB population (Nava et al." 
[2006 ). As sample sizes slowly increase, addressing these issues re- 
mains critical for the correct inference of intrinsic GRB properties 
and thus their use as cosmological probes. 

In this paper, we concentrate on a long-standing empirical 
relation that initially suggested a possible Cepheid-like correla- 
tion between gamma-ray variability and peak luminosity of GRBs 
(Reichart et al. 2001; Fenimore & Ramirez-Ruiz 2000). 

iReichartetaljiioOlh (hereafter ROl), using a robust measure 
V of the GRB variability, for a sample of 13 GRBs with known red- 
shift, found that in the GRB rest frame this measure is correlated 
with the GRB peak luminosity L. In the log L-log V plane the cor- 
relation was modelled with a linear function log L = m log V + b 
with the slope of the straight line m = S.S^J'g and a sample 
variance along V of aiogv = 0.18, both parameters being ob- 
tained with the method described by Reichart 1 2001) (hereafter 
called Reichart method). This method was proposed to fit data sets 
affected by a sample variance in addition to the statistical vari- 
ance (called "intrinsi c variance") of each data poi nt. Recent l y, firs t 
iGuidorzi et ail 12005) (hereafter GFM05) and then Guidorzi 1 2005) 
(hereafter G05) tested this correlation first using an extended sam- 
ple of 32 GRBs with known redshif t (GFM05), and the n with 551 
GRBs detected by CGRO/BATSE iPaciesas etal."l999|) for which 
a pseudo-redshift was derived by G 05 exploiting the s pectral lag- 
luminosity correlation (Norris et al. 2000; Norris 2002; B and et alJ 



12004) . In both cases, the correlation was confirmed and in the 
log V-log L plane the slope of the straight line was found much 
lower than that derived by ROl (m = 1.31q;4 derived by GFM05; 
m = 0.85 ± 0.02 derived by G05). It was also found that, with 
the sample variance neglected, the straight line did not provide a 
good description of the data (x^/dof — 1167/30 and x^/dof = 
4238/549 for the samples considered by GFM05 and G05, respec- 
tively). Neglecting the sample variance was correctly questioned 
by^eichart & Nysewander ( 2005) (hereafter RN05), who however 
show that, also using the extended data set of GFM05, the m slope 
is given by m = 3.4+Qg, still in perfect agreement with the origi- 
nal value found by ROl and in strong disagreement with the value 
found by GFM05. They attribute this disagreement to the fact that 
GFM05 do not include among the parameters of the fit the sam- 
ple variance, which they estimate to be aia^v = 0.20 ± 0.04. In 
an unrefereed note, lReichart| i2005h takes issue with the most re- 
cent paper by G05, who confirms the results previously found by 
GFM05. 

In this paper, we discuss the Reichart method and compare the 
results with those obtained following the treatment by D'Ago stinL 
(2005) (hereafter called "the D'Agostini method"), which deals 
with the same problem of fitting data points affected by extrinsic 
scatter in addition the intrinsic uncertainties along both axes. We 
show with numerical simulations that the Reichart method provides 
an inconsistent estimate of the m slope in the specific case of the 
log V-\og L data set, we discuss the likely reason for this incon- 
sistency. We show that our original results are substantially con- 
firmed by the D'Agostini method even taking into account the sam- 
ple variance. Finally we discuss the subsequent implications for the 
inferred physics of GRB central engines and their relativistic out- 
flows. The usefulness of the so-called variability/peak luminosity 
correlation is discussed in the broader context of recently discov- 
ered correlations between other observed and derived properties of 
GRBs. 



2 THE REICHART METHOD 

The Reichart method has as starting point the well known Bayes 
theorem that, for inference of physical parameters, is widely 
discussed in several text books (see, e.g., D'Agostini 2003). It 
states that the conditional probability that a set of parameters 6= 
0\, 02, ■■■On takes a particular value 0", given a data set D, whose 
values depend on 6 (for instance, N independent observations of a 
quantity X), is given by: 



p{9\DI) = 



p{e\i) v{D\ei) 

Pirn 



(1) 



where / is the available prior information, p{d\I) is the probability 
distribution of on the basis of the information /, and p{D\0I) is 
the conditional probability of getting the measured data set given 
the value 9" and the information /. The probability distribution 
p{0\I) is called prior probability, while the probability distribu- 
tion p(_D|0J) is called likelihood function. The probability p{D\I) 
is introduced as normalization factor. In the case that no prior in- 
formation is known, p(0\I) is generally assumed to be a uniform 
distribution and the range of the possible values of 6 are those log- 
ically allowed. In this case the posterior probability and the likeli- 
hood function are equivalent. 

Reichart (2001) concentrates on the derivation of the prior 
probability p{0\I) in the special case in which for two of the pa- 
rameters, X and y (in our specific case x = log V and y = log L), 
we have a set of pairs of measurements, from which it appears 
that the two are correlated with y — yc{x;6m), were 6m, with 
m = 1,2,..., M, are M intermediate parameters that describe the 
curve yc{x; 6m) (in our case y — mx + b, M = 2). Correctly 
the scatter of the parameter values x and y around the curve is as- 
sumed to be due partly to the measurement errors (intrinsic scat- 
ter) and partly due to weaker dependences of either parameters x 
or y on other, yet unmeasured, and even yet unknown, parameters 
(extrinsic scatter or sample scatter). Both the scatters for x and y, 
are assumed by Reichart to be normally distributed and uncorre- 
lated, with unknown standard deviations and ay of the extrinsic 
scatter of x and y, respectively. The conditional probability of the 
values of the parameters dm,<yx, Oy given the measured set of data 
points and their uncertainties is thus derived under simplifying as- 
sumptions (see Section 2.2.2 of Reichart 2001), among which that 
the curve y = yc(x; dm) can be approximated by a straight line 
(y ~ yt,n + St.n{x — xt.n)). The rcsult is the following (eq. 43 of 
Reichart 2001): 



p{dm:Ox,CJy\Xn,yn-,Ox,n,Oy,„) ~ V 1 + ^t.n Cn \yu,yt,n + 



+ St,„(a;„ - Xt.n), \J cr.y + a|,„ + + a|,„)] (2) 

where cr^ „ and (j'y^„ are the intrinsic variances of the A'^ pairs 
of data points Xn and yn (n = 1,2, TV), respectively; 
{xt.n, yt.n) is the point on the curve y — ydx; 6m) which max- 
imises the two-dimensional Gaussian Gn{x,x„, -^/crl + cr|,n) x 
Gn{y,yn, \J<yl^ o-g,„). 

In order to simplify the derivation of the prior probability in 
this special case, eq.|2|is assumed to be a likelihood function and 
the maximum likelihood method is applied to constrain the values 
of the intermediate parameters 6m, and Oy and the uncertainty 
in the values of 6m ■ 
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Figure 1. Variability vs. Peak Luminosity for tlie 32 GRBs with known 
redshift considered by GFM05. Also shown are the best fit curves and Itr 
regions. Solid lines: best fit results with the Reichart method. Dashed lines: 
results obtained by GFM05. Dashed-dotted lines: best fit results with the 
D'Agostini method. 
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Figure 2. Variability vs. Peak Luminosity for the 551 BATSE GRBs with 
pseudo-redshift considered by G05. Also shown are the best fit curves and 
l(j regions. Solid lines: best fit results with the Reichart method. Dashed 
lines: results obtained by G05. Dashed-dotted lines: best fit results with the 
D'Agostini method. 



2.1 Application of the Reichart method to the 
Luminosity- Variability correlation 

We now apply the Reichart method in the specific case of the 
variability-luminosity correlation test, showing that we are capa- 
ble of reproducing the results obtained by ROl and RN05. In the 
case of this test, N pairs of measured values of Vn and Ln (one 
per each burst) are available, and, in the eq.|2| x„ — log {V„/V), 
Un ~ log (Ln/L) (V and L being the correspondent median val- 
ues), ax,n = o"iogv„ and ay.n = criogL,,- Having modelled the 
correlation by a straight line: 

y{x) — mx + q (3) 
from eq.|2| it is possible to derive the log-likelihood function: 

logp(m, g, (Jx,ay\{xi,yi,(Jx,i,cjy^i}) = 

N 



+ 



ay +m? a'l - + m? 



(4) 



in which, without loss of generality, it is possible to assume either 
ax ~ 0, or, alternatively, ay = (both appear only in the term 

2 I 2 2n 

ay+m a J. 

Assuming a^ = in eq. |4| with the maximum likelihood 
method the free parameters m, q, and ay can be derived. Alter- 
natively, if one assumes ay — 0, the value for ax can then be de- 
termined from ax = ay/m, where ay is the value obtained in the 
previous case under the assumption of null ax ■ 

For our extended sample of 32 GRBs with known redshift 
(GFM05), we find that at 90% confidence level, the best-fitting pa- 



rameter values are the following: m 

1+0.77 



3.8 



+2.8 



0.07 



+0.32 



aiogL 



0-93Io:29 (alternatively, ax = criogv ^ 0.24) 



(solid lines in Fig.0. As can be seen, all these values are in ex- 
cellent agreement with those derived by RN05 (they report the 
la uncertainty). Similarly, for the sample of 551 BATSE GRBs 



with pseudo-redshift (G05), we find m = 3.5to'4> Q 
ay = 1.2X1:°:", or ~ 0-35 (solid lines in Fig.|2j. All these 
results are also summarised in TableQ 



0.15 



+0.9 



3 AN UNBIASED METHOD 

iD'Agostini i2005t) addressed the same problem described in the 
previous section, i.e. how to perform a linear fit between two data 
sets with errors on both axes and with an extra variance. Simi- 
larly to the Reichart method, the D'Agostini method is based on 
the parametric inference typical of the Ba yesian approach. H ow- 
ever, the resulting log-likelihood derived bv lD'Agostinilj200^ (see 
eq. 35 and 52 therein) differs from that by Reichart, reported in|2| 
for just one term: in the D'Agostini likelihood function the term 
(1 -f m^) is just 1. For a detailed description of the D'Agostini 
method and of its derivation we address the reader to the original 
paper. We will demonstrate that, unlike the Reichart method, the 
likelihood function derived by D'Agostini provides unbiased esti- 
mates of the unknown parameters. 

According to D'Agostini (2005), the log-likelihood function 
for the case here considered, is considered is given by 

logp(m, q, ax, ay\{xi,yi, ax,i,ay,i}) = 



N 

E 

i = l 



[^°^(2^(a2 + m2 



a'i + a'^ - + m? 



if 



+ m- 



+ a^ - + m2 



(5) 



Using this equation, in the case of the 32 GRBs of GFM05, at 90% 
confidence level we find the results reported in TableQ In particular 
we find a slope m = 1.7 ± 0.4 against a value m = 3.81^ * found 
with the Reichart method. Likewise, using the sample of 551 GRBs 
with pseudo-redshift of G05, we find the best fit parameter values 
reported in TableQ m = O.SSlg iJg. The best-fitting power-laws 
obtained with the D'Agostini method and the Reichart method are 
compared in Fig.sQand|2|for the two data sets, respectively. 

In Table Q we report the also the best-fitting parameters ob- 
tained using the Least square fit in the case of data in the case of 
data affected by errors on both axes ("f itexy" tool. IPress et alJ 
( 1993)), but not with no extra variance. We report these results to 
examine how the best-fitting parameters are affected when the sam- 
ple variance is taken into account. 

As can be seen from Table \\\ the Reichart method and 
D'Agostini method give different best fit parameter values, espe- 
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Table 1. Best-fitting parameters obtained with different methods for the GRB samples used by GFM05 and G05. The confidence intervals are at 90%. 



ORB Set Method m q cry ctx 

(o-^ = 0) (cry = 0) 



32 (from GFM05) 


Reichart 


3.8+^-^ 


'^•'^'-0.33 


93+0-77 

u.»o_Q 29 


~ 0.24 


32 (from GFM05) 


D'Agostini 


^ 7+0:4 

'-■ ' -OA 


-0.19 


0.58«:1^ 


~ 0.34 


32 (from GFM05) 


f itexy 


1.9±0.l(=') 


0.14±0.02(=') 






551 (fromG05) 


Reichart 


q c+o.e 


U.iO_Q g 


1 91+0.19 
^•^^-0.15 


~ 0.35 


551 (fromG05) 


D'Agostini 


0.881°;- 


01+0 03 

^^•"^-0.03 


n 6"s+o-o-* 


~ 0.74 


551 (fromG05) 


f itexy 


1.37±0.02(=') 


0.09±0.0l(^) 







1(7 confidence interval. 



cially the value of the slope m. The Reichart method yields steeper 
slopes even that those found with the D'Agostini method, whose 
results are consistent with those originally obtained by GFM05 and 
G05 (m = l.StJS j, GFM05; m = 0.85 ± 0.02, G05), even though 
the fit, in the latter case, was found unac ceptable (ver y high Xr)- 
Also with the 'f itexy' algorithm (Pre ss et alJll993h . the slopes 
obtained are much shallower than the correspondent obtained with 
the Reichart method. 



4 COMPARATIVE ASSESSMENT OF ANALYSIS 
METHODS 

In order to test the consistency of the estimate of the m param- 
eter obtained with each method, we performed various numerical 
simulations. 

Assuming the existence of a linear correlation between y = 
log {L/L) and x — log (V/V) with prefixed parameters (m, q, 
ax, o-y), using random number generators, proper data sets were 
produced and the data analysis was performed following both the 
Reichart and the D'Agostini methods. The best-fit parameters were 
thus compared with the corresponding input values, also called 
"true values". We assumed ay = without loss of generality, as 
discussed in Sec. 12. II We call ax,trua the value assumed to gener- 
ate the populations of random {xi,yi) points, in order to distinguish 
it from the best-fitting values derived with each method. Likewise, 
we denote mtruo the correct value for m. 

From our tests it resulted that there are two different regimes, 
mainly depending on the value of the ratio between the extrinsic 
variance a^ and the total variance of the sample along x, a^ ^. 

Concerning the distribution of intrinsic uncertainties along x 



and y, ax,i and ay^i (i — 1, 
two different cases: 



, A''), respectively, we considered 



(i) the {ax,i} have been generated from the distribution of 
the intrinsic uncertainties along x of the sample of 32 GRBs of 
GFM05. Figure|3| shows the two populations of ax.i as a function 
of X in two cases: the real case of the 32 GRBs of GFM05 and a 
random population of 200 points. The same was done along y to 
generate the distribution of cry, i; the result is shown in Fig.|4| 

(ii) the same as the previous case, but with a population of in- 
trinsic uncertainties {ax,i} a factor 1/10 times as big as those gen- 
erated following the previous method. The reason is to study the 
regime - << a^^f 

For each simulated data set, we generated 1000 points accord- 
ing to a uniform distribution in x in the range — 2 < x < 2. We 
then selected 200 points in the range —0.5 < a:: < 0.5. This choice 
is because we want a range along a; of 1, which is the same as that 
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Figure 3. Intrinsic uncertainties a^j asa function of x for two populations: 
200 random simulated points (crosses) and the 32 real GRBs of GFM05 
(filled circles). 
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Figure 4. Intrinsic uncertainties cry,,; as a function of y for two populations: 
200 random simulated points (crosses) and the 32 real GRBs of GFM05 
(filled circles). 



of the 32 GRBs of GFM05, ranging in -1.6 < log 1^ < -0.6 
(Fig.D. 

Concerning o-j^^truc, we studied the dependence of m and a^ 
as a function of ax,tiuc in the range < o"a:,truc < 0.5. 

• Case a. We assumed mtruc ~ 1-7, gtruc = (similar to that 
derived with the D'Agostini method for the GFM05 sample), and 
a distribution of intrinsic uncertainties along x and y, {a^.i} and 
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Figure 5. m as a function of era;, true. Solid line show.s mtruc = 1-7 (Case 
a). Empty squares con'espond to the D'Agostini, tilled circles to the Re- 
ichart method. Error bars show the 90% contidence interval. 
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Figure 6. ax estimated with either method as a function of the true value 

0-a;,true (CaSe a). 



{oy^i}, given by item (i) (Figs.|3|and|4j. The estimated m as a func- 
tion of the true extrinsic scatter along x, a^.tiuc, is shown in Fig.|5] 
for either method. Apparently, when the ratio a-x.tiuc/o'x,t > 0.4, 
m estimated with the Reichart method is significantly overesti- 
mated, while that derived with the D'Agostini method turns out 
to be a more consistent estimator. 

Figure |6| shows ax derived with each method as a function of 
the true value CTa:,truc. In this example, for (Jx.tTuc/o'x,t > 0.4, the 
Reichart method seems to underestimate the extrinsic scatter ax 
(Fig-EJ and overestimate the slope m (Fig.|5j; these two aspects 
appear to be correlated. FigureQshows the case of ax.tmc = 0.34, 
similarly to the result of the D'Agostini method applied to the 
GFM05 sample. In that case, the estimated slopes result: m = 
3.6lg 4 (Reichart), m = 1.7 ± 0.2 (D'Agostini). The total scat- 
ter along X is ax,t = 0.285. 

• Case b. We adopted the same values for mtruc and gtruc as 
those of Case a, but with intrinsic errors 10 times smaller along 
both axes, as described by (ii). The results are basically the same 
as Case a, confirming that when ax, true becomes comparable with 
ax,t the Reichart method overestimates m and underestimates ax- 
When (7a;, true = 0.34 it rcsults: m = 3.7 ± 0.5 and ax ~ 0.22 
(Reichart), m = 1.6 ± 0.2 and ax ~ 0.35 (D'Agostini). 



Figure 7. Sample of 200 simulated points generated assuming '^truc — 
1.7, QtTuc — (Case a) and cr^;, true — 0.34. Dashed lines show the best- 
fit and 1(T region obtained with the Reichart, while dotted hues correspond 
to the D'Agostini method. Solid thick lines show the true power law and 
dispersion. 
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Figure 8. m as a function of (Ta;,true. Solid line shows mtruo = 3.8 (Case 
c). Empty squares correspond to the D'Agostini, filled circles to the Re- 
ichart method. EiTor bars show the 90% confidence interval. 



• Case c. We assumed mtrue = 3.8, qtrue ~ (similar to that 
derived with the Reichart method for the GFM05 sample), and a 
distribution of intrinsic uncertainties along x and y, {ax,i} and 
{ay.i}, given by item (i) (Figs.|3|and|4}. The estimated m as a func- 
tion of the true extrinsic scatter along x, cri;,truc, is shown in Fig.|8| 
for either method. Similarly to what found in Case a, when the ra- 
tio o-a;,true/o"a:,t > 0.2, 771 estimated with the Reichart method is 
significantly overestimated, while that derived with the D'Agostini 
method proves to be a more robust estimator. 

Figure |9| shows the case of ax,tTuc ~ 0.24, similarly to the re- 
sult of the Reichart method applied to the GFM05 sample. In that 
case, the estimated slopes result: m = and ax ^ 0.17 

(Reichart), m = S-YILo ^ and ax ~ 0.22 (D'Agostini). The total 
scatter along x is ax.t = 0.286. 

• Case d. We adopted the same values for mtrue and gtruo as 
those of Case c, but with intrinsic errors 10 times smaller along 
both axes, as described by (ii). The results are consistent with 
those of Case c, confirming that when CTa:,true becomes compara- 
ble with ax,t the Reichart methods overestimates m and underes- 
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Figure 9. Sample of 200 simulated points generated assuming TTitruc — 
3.8, (jtrue = (Case c) and = 0.24. Dashed lines show the best-fit 
and Icr region obtained with the Reichart, while dotted lines con'espond 
to the D'Agostini method. Solid thick lines show the true power law and 
dispersion. 



timates cfx- When (7a;,truc ~ 0.24 it results: m = 6.4 ± 0.6 and 
(Jx ~ 0.18(Reichart), m = S.OtJJlg and ~ 0.22 (D'Agostini). 

To better appreciate what happens when the total scatter Ox,t 
of the sample becomes comparable with the extrinsic scatter Ox, as 
it might be the case when truncation effects are at work, in Fig. 1101 
we show the the entire population of 1000 random points generated 
in Case a, when (Tx,true = 0.34, in the range —2 < x <2 and the 
selected subsample of 200 points in the range —0.5 < x < 0.5. 
The total scatter ax,t is 1.39 {ox,true/ <Jx,t = 0.24) for the entire 
sample of 1000 points and is 0.285 ((Tx,truE/ o-x,t = 0.84) for the 
subsample of 200 points. 



5 DISCUSSION 

The results obtained from the numerical simulations show that, 
when the sample variance cr^ truo is significantly lower than the to- 
tal variance o"^^t, both methods provide a consistent estimate of the 
m slope. Differently, when cr^^truc is comparable with cr^ t, the Re- 
ichart method overestimates m, while the D'Agostini method still 
gives an estimate of m compatible with the correct value mtiuc. 

The simulations have been performed adopting the same dis- 
tributions of intrinsic errors along both axes, {(yx,i} and {(7y,i\ 
(j = 1, . . . , iV), as those of the sample of 32 GRBs of GFM05. For 
each simulated data set, we selected 200 points within the range 
—0.5 < X < 0.5, similarly to the range covered in log [V/V) by 
the 32 GRBs of GFM05. Typical values of the total scatter Ox,t for 
the simulated data sets are 0.28-0.29, to be compared with that of 
the 32 GRBs of GFM05, which is 0.32. The results obtained are 
confirmed even when we adopt the same extrinsic scatter and 10 
times smaller intrinsic errors along either axis. The most important 
factor is the ratio between the extrinsic scatter, (Jx,trvLc, and the total 
scatter, ax,t- 

The two methods differ when the range along x is compara- 
ble with the extrinsic scatter along the same axis. This is shown 
in Fig. ^1 the full population clearly follows the true slope of 
W2true = 1-7 shown by the cyan line. However the two methods, 
applied to the 200 points selected in the range —0.5 < x < 0.5, 
for which the ratio (Tx,trvLc/ crx,t = 0.84, give different results, the 



Reichart method (green) overestimating m, while the D'Agostini 
method (yellow) is still consistent with the true slope. We do not 
know yet whether the V — L correlation presently known ex- 
tends along V over a wider range, the samples considered by 
Reic hart et alJ ^001) and GFM05 being the result of truncation 
effects, or it does not. Either way, however, the total variance is 
comparable with the extrinsic variance of the correlation and in 
this case, the D'Agostini method turns out to be a more reliable 
estimator of m than the Reichart. 

Similar results to those presented here and by GFM05 and 
G05, i.e. values of m around 1-2, have already been mentioned 
by Llovd-Ronning et al. 1 200 2): in fact, therein we read that a 
more recent analysis by iFenimore & Ramirez-Ruij llOOO) mak- 
ing use of a different definition of variability, based on a differ- 
ent degree of smoothing of the light curve, led to an estimate of 
m = 1.57± ~ 0.5. This supports the view of a not well estab- 
lished value of power-l aw index of the V — L c orrelation, affected 
by a significant scatter iLlovd-Ro nnin g et al. 2002). 

In another recent paper bv fLi & Paczvnski ^2006) a new def- 
inition of variability is considered, based on a different smoothing 
filter of the light curves. In this work, they find a slope of ~3, ap- 
plying the f itexy procedure, and no extra variance is needed, as 
the correlation appear s to be tight enough. These results as well as 
the recent analysis bv lFenimore & Ramirez-Rui3 i2000l) point to a 
strong dependence of m on the kind of filter adopted to obtain a 
smoothed version of the light curve, with respect to which variabil- 
ity of the original light curve is computed. 

As far as th e definition of variability introduced by 
iReichart et alj i200lh is regarded, we have investigated reasons for 
the overestimation of m given by the Reichart method and found 
that this is likely an incorrect derivation of the likelihood function 
by Reichart (2001) (eq.|5J- This function was derived by integrating 
the weighted product of Gaussian density functions over, among 
other variables, the element ds'^ — {dx'^ + dy^) (eq. 28 of Re- 
ichart 2001), where y — yc{x,9m) is the equation which describes 
the relation between the parameters x and y (see Section|2}. The 
integration over ds^ = {dx^ + dy^) is not correct: dy is homoge- 
neous with m dx and not with just dx. The integration on ds in- 
troduces the sum 1 + in the numerator of the first term at right 
hand side of eq.|4| which is the only difference between the likeli- 
hood functions of the two methods. The inconsistency of the term 
ds mentioned above is transmitted to the factor 1 + m^: m can- 
not be added to 1 tout court for the same reason explained above 
and discussed bv lD'Aeostin i' ('2005'). The effect of this factor is the 
overestimation of m. This is shown by the fact that, by deleting 
m^, the estimate of m becomes consistent and the res ulting likeli- 
hood function becomes coincident with that derived bv lD'Agostinil 
(2005) for the same statistical inference problem, reported here in 
eq.|5l 

The procedure adopted by GFM05 and G05 as well as the 
usage of the f itexy procedure reported in TableQdoes not take 
into account the sample or extrinsic variance. As a result they find a 
high assuming a straight line for the correlation between log L 
and log V. However, if, in spite of that, a straight line is assumed, 
the line slope m value found is very close to the true one. 

We conclude that the results reported by ROl and RN05, who 
applied the Reichart method, give an overestimated index m of the 
power law (L oc V™") function used to describe the correlation be- 
tween peak luminosity and variability of the GRBs included in the 
sample of ROl, GFM05 and G05. We ascribe the different results 
obtained with the D'Agostini and Reichart methods, when applied 
to the samples of 32 GRBs of GFM05 and 551 GRBs of G05, to 
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Figure 10. Full population of 1000 simulated points (crosses) from Case a, with mtruc = 1-7 and Caj^truc 
subsample of 200 points —0.5 < x < 0.5, both fitting methods have been applied to in Case a. 
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the comparable total variance cr^. ^ and extrinsic variance a^. The 
D'Agostini method appears to provide the most correct estimations 
of the slope and extrinsic scatter of either sample. The best-fitting 
power laws derived by GFM05 and G05 are more consistent with 
the results of the D'Agostini than with the Reichart method. 

The implications of different power-law index and extrinsic 
scatter describing the V ~ L correlation are twofold: the physical 
interpretation and its usage as a luminosity estimator. 



5.1 Implications on the interpretation 

Variability in the GRB light curves is thought to be produced 
mainly above the photospheric radius at which the relativistic out- 
flow becomes optically thin to electron scattering. In addition to 
this photosphere, there would be a second photospheric radius be- 
yond which shocks are optically thin to pairs too and highly vari- 
able light curves originate. Meszaros et al. (2002) interpreted the 
V — L correlation as due to the effects of this pair-forming pho- 
tosphere and the jet opening and/or viewing angles, when breaks 
are observed in the optical consistent with the jet interpretation. 
The power-law m derived by those authors is m = 6p/{5q — 2p), 
where they assumed L oc 9~^, p ~ 2 — 2.5 (6 can be either the 
jet opening angle or the viewing angle for a universal jet pattern), 
Fm oc 9~'' (T,n is the minimum value of the Lorentz factor distri- 
bution of the shells). Following Meszaros et al. (2002), if one takes 
p — q, the expected power-law index is m = 2, which is consis- 



tent wit hin 90% CL with the va lue we derived with the D'Agostini 
method. iMeszaros et alj i2002l) assume d p = 5/2 and q = 2 to 
ex plain the value originally found b y iReichart et alj j2 00ll) and 
by Fe nimore & Ramir ez-Ruiz (2000) of m ~ s iKobavashi et all 
(2002) started from similar assumptions: F cx 0"', i.e. a wider jet 
has a lower Lorentz factor, equal mass colliding shells, a uniform 
D/c distribution in logarithmic spaces in the interval 1 ms-1 s (D 
is the separation betw een shells). Using a simplified version of the 
IReichart et alJ i200lh variability they tried to reproduce the V — L 
correlation with numerical simulations and found that q = 1 seems 
to better account for m values in the range of 3—4 , as fo und by 
IReichart et ail <200ll) and lFenimore & Ramirez-Ruidiioool) . How- 
ever, from the same numerical results it turns out that smaller val- 
ues of m, compatible with the results obtained with the D'Agostini 
m ethod, are given by ass uming q = 2. As already pointed out 
bv lKobavashi et all J2002l). this is more consistent with the value 
g = 8/3 derived by Salmonson & Galama ('2002*), based on the 
study of the anticorrelation between peak luminosity and jet break 
time. From the results of lKobavashi et alj i2002t) . we infer that the 
results of the D'Agostini best-fitting parameter m imply a stronger 
dependence of the Lorentz factor on the opening angle as well as a 
smaller Lorentz factor normalisation. 



Likewise. lRamirez-Ruiz & Llovd -Ronning ( 2002) studied the 
effects of varying the energy per unit solid angle in the collimated 
outflow on the correlatio ns between variability, peak luminosit y 
and spectral peak energy ^Llovd-Ronning & Ramirez-Ruizll20o3) . 
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They found that these correlations are accounted for if luminos- 
ity per unit solid angle L{9) (where 6 can be either the view- 
ing angle for a structured jet pattern or the jet opening angle), 
strongly depends on the Lorentz factor T{6). while the effect of 
changes in the baryon loading of the wind is less relevant, as it 
does not account for the range of variability observed. In their 
figure 3,|piamirez-Ruiz & Llovd-Ronnins (2002) show the results 
in the V ~ L space of some numerical simulations for different 
bulk Lorentz factors: values for m in the range 1.5-2, i.e. like 
those found by GFM05 and those derived in this work with the 
D'Agostini method, are clearly compatible with the results of these 
simulations. 

The smaller value of m we found with the D'Agostini method 
than with the Reichart constrains differently the parameters m de- 
pends on. Moreover, the bigger extrinsic scatter of oio^v, 0.34 with 
the D'Agostini, 0.24 with the Reichart method, means a bigger 
scatter in the distribution of some of these parameters. 

5.2 Implications on its usage as luminosity estimator 

After the discovery o f the V — L correlat i on b y 
iFenimore & Ramirez-Rui3 <2000t) and iReichart et all <200ll) . 

a numb er of papers used it as a luminos ity estimator for several 
studies. IFenimore & Rami rez-Ruiz' teoool) applied it to derive the 
luminosity distribution 0(-L) and the G RB form ation rate as a 
function of redshift oiz^. iLlovd-Ron ning et al. (2002) studied the 
bivariate distribution of luminosity and redshift ^{L,z) using the 
pseudo-redshifts obtained from a parametrisation of the V ~ L 
correlation with m — 1.57± ~ 0.5. They found a correlation be- 
tween luminosity and redshift parametrised as L oc (1 + z)^'^*'^ *. 
Different valu es of m from ^ 1 to ^ 3 d o not change their results 
qual itatively ^lovd-Ro nning et all2002l) . 

iLazzatil 2002) exploited the V — L correlation to study the 
rest-frame power spectral density (PSD) in time domain for hun- 
dreds of BATSE bursts. In particular, the goal was to study the ef- 
fect of photon scattering on variability. A correlation between the 
break frequency in the rest-frame PSD and variability was found, 
based on m ~ 3. A smaller value of m would mean that a GRB 
with a given variability, whose dependence on z is small, would be 
less luminous, i.e. nearer, thus with a corresponding smaller break 
frequency in the rest-frame PSD. This relation would have to be 
rescaled consequently. 

More recently, other tighter correlations became popular, such 
as between the rest-frame peak energy Ep and the isotropic- 
equivalent gamma -ray released 

jAmatietalJl2002h . between Ep and the coUi mation-corrected 
SMnma-ray released energy E-^, Ep oc E^''^ iOhirlanda et alJ 
|200J). Another relation between Ep, £^,150 and the rest-frame 
break time % of the optical afterglow light curve has been dis- 
covered recently iLiang & Zhan3l2005h . which requires fewer as- 
sumptions than those of the Ghirlanda relationship. The same re- 
lation as the Amati holds with the same index of ~0.5 when 
E-fjso is replaced by the peak luminosity L-^aso, or simply L 
JVonetoku et al . 2004; G hirlanda et alJboOSi) . Recently, the Amati 
relation has been confirmed for a sample of 53 long GRBs and a 
fit with the D'Agostini method gives an extrinsic logarithmic s cat- 
ter of 0.15 over 3 orders of magnitude over Ep ( Amati '200f?). In 
this case both methods give compatible results and this is consis- 
tent with the fact that ax « cTx.t- Summing up, in the case of the 

Amati relation, truncation effects are negligible^ 

Interestingly, iLlovd -Ronning & Ramirez -Ruid i2002h stud- 
ied the correlation between V and Ep before the discovery 



of the Amati relationship. Assuming different alternative val- 
ues for m, ranging from 2.2 to 5.8, as originally found by 
IFenimore & Ramirez-RuizI i2000() . they found a correlation, Ep oc 
V\ with 5 ranging from 0.4 to 1.15 . The Amati relation, whose 
extrinsic scatter is small with respect to that of the V — L corre- 
lation, is more reliable as a luminosity estimator. Hence, from the 
Amati relation we should expect m = (5/0.5, so in the range 0.8- 
2.3, which is consistent with the results of the D'Agostini method 
as well as those of GFM05 and G05. 



6 CONCLUSIONS 

We applied both Reichart and D'Agostini methods to the sam- 
ples of 32 GRBs with known redshift and 551 BATSE GRBs with 
pseudo-redshift considered by GFM05 and G05, respectively. The 
goal was to estimate the slope m as well as the scatter of the power 
law describing the correlatio n between variability and peak lumi- 
nosity originally presented bv lReichartetalJj200l lh. Both methods 
account for an extra variance in addition to the intrinsic affecting 
each single point. From simulations, we found that when the sam- 
ple variance is comparable with the total scatter along the same 
axis, (T^ t, the Reichart method tends to overestimate m, while the 
D'Agostini still estimates it correctly. When the sample variance is 
negligible with respect to the total variance, the two methods give 
consistent results. In the specific case of the V — L correlation, 
the two variances are comparable: in the case of the 32 GRBs of 
GFM05, it is ax,t = 0.32, = 0.24 (Reichart) or = 0.34 
(D'Agostini). This explains the discrepancy between the two meth- 
ods. We showed that the D'Agostini method is a reliable estima- 
tor of m in this regime, whereas the Reichart is no more. In par- 
ticular, the best-fitting value for m obtained with the D'Agostini 
method is 1.7 ± 0.4 and 0.88lo;i3 at 90% confidence level for the 
32 GRBs of GFM05 and for the 551 GRBs of G05, respectively. 
These values are significantly smaller than those obtained with the 
Reichart method, which are consistent with previous results of ROl 
and RN05. These results hold as far as the definition of variability 
given by ROl is assumed. 

Alternatively, from other definitions of variability based on 
different kinds of filter used to smooth the light curves, it seems 
to be possible to obtain a range of values for m from '^l to ^^3 
jFenimore & Ramirez-Ruizl200(llLi & Paczvnskj2006h . 

One of the possible implications of a smaller value of m than 
originally found by ROl on the interpretation of the V — L correla- 
tion is that, in the jet emission scenario, we would expect a stronger 
dependence of the Lorentz factor of th e expanding shells on the jet 
openmg ang le teobavashi et alj2002h . 

Finally, other more recent and tighter correlations, such as 
the Amati, Ghirlanda, Liang & Zhang relations appear to be 
more reliable luminosity estimators than the V ~ L one. In 
particular, our results of values of m in the interval 1-2 ob- 
tained with the D'Agostini method, appear to be consistent with 

two independent relations: Ep ex. (0.4 < 5 < 1.15) 

( iLlovd-Ronning & Ramirez-Ruizll2002h . and the equivalent ver- 
sion of the Amati relation with the isotropic-equivalent peak lu- 
minosity instead of the isotropic-equivalent total released energy, 
Ep oc L"'^. By combining the two, one expects m — 5/0.5, in 
agreement with the results reported here as well as those reported 
by GFM05 and G05. 

Finally, the increasing number of GRBs with spectroscopic 
redshift detected by Swift ( Gehrels et al. 2004) will help extend the 
range in V and better constrain the power-law fit of the V — L cor- 
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relation, with the benefit of a homogeneous data set of fight curves 
afi detected with the Burst Alert Telescope (BAT) onboard Swift. A 
thorough test of the V ~ L correlation with BAT data is in progress 
(Rizzuto et al., in preparation). 



ACKNOWLEDGMENTS 

We thank the anonymous referee for useful comments. C.G. thanks 
S. Kobayashi for useful discussions. C.G. acknowledges his Marie 
Curie Fellowship from the European Commission. CO., F.F., 
E.M., F.R. and L.A. acknowledge support from the Italian Space 
Agency and Ministry of University and Scientific Research of Italy 
(PRIN 2003 on GRBs). A.G. acknowledges the receipt of Marie 
Curie European Re-integration Grant. C.G.M. acknowledges fi- 
nancial support from the Royal Society. This research has also 
made use of data obtained from the HETE2 science team via 
the website http://space.mit.edu/HETE/Bursts/Data and BATSE, 
Konus/WWD and BAT/Swift data obtained from the High-Energy 
Astrophysics Science Archive Research Center (HEASARC), pro- 
vided by NASA Goddard Space Flight Center 



Reichart D.E., Lamb D.Q., Fenimore E.E., Ramirez-Ruiz E., 

Cline T.L., Hurley K., 2001, ApJ, 552, 57 (ROl) 
Reichart D.E., 2001, ApJ, 553, 235 (RM) 

Reichart D.E., Nysewander M.C., 2005, ApJ, submitted 

(astro-ph/0508111 1 (RN05) 
Reichart D.E., 2005, astro-ph/0508529 
Salmonson J.D. & Galama T.J., 2002, ApJ, 569, 682 
Yonetoku D., Murakami T, Nakamura T, Yamazaki R., Inoue 

A.K. & loka K., 2004, ApJ, 609, 935 

This paper has been typeset from a TjiX/ E'lgX file prepared by the 
author. 



REFERENCES 

Amati L. et al., 2002, A&A, 390, 81 
Amati L., 2006, MNRAS, submitted I astro-ph/0601553'1 
Band D.L., Morris J.R, Bonnell J.T, 2004, ApJ, 613, 484 
D'Agostini G., 2003, Bayesan reasoning in data analysis: A crit- 
ical Introduction, World Scientific Publishing. 

D'Agostini G., 2005, physics/051 1 182 

Fenimore E.E., Ramirez-Ruiz E., 2000, astro-ph/00 04T76| 
Gehrels N. et al., 2004, ApJ, 611, 1005 

Ghirlanda G., Ghisellini G. & Lazzati D., 2004, ApJ, 616, 331 
Ghirlanda G., Ghisellini G., Firmani C, Celotti A. & Bosnjak Z., 

2005, MNRAS, 360, L45 

Guidorzi C, Frontera F., Montanari E., Rossi F., Amati L., Gom- 
boc A., Hurley K., Mundell C.G., 2005, MNRAS, 363, 315 
(GFM05) 

Guidorzi C, 2005, MNRAS, 364, 163 (G05) 

Haislip J.B. et al., 2006, Nature, 440, 181 

Kawai N. et al., 2006, Nature, 440, 184 

Kobayashi S., Ryde F & MacFayden A., 2002, ApJ, 577, 302 

loka, K. & Nakamura, T, 2001, ApJ, 554, L163 

Lazzati D., 2002, MNRAS, 337, 1426 

Li L. & Paczynski B., 2006, MNRAS, 366, 219 

Liang E. & Zhang B., 2005, ApJ, 633, 61 1 

Lloyd-Ronning N.M., Fryer C.L. & Ramirez-Ruiz E., 2002, ApJ, 
574, 554 

Lloyd-Ronning N.M. & Ramirez-Ruiz E., 2002, ApJ, 576, 101 
Meszaros P., Ramirez-Ruiz E., Rees M.J., & Zhang B., 2002, ApJ, 
578, 812 

Nava L., Ghisellini G., Ghirlanda G., Tavecchio F. & Firmani C, 

2006, A&A, 450, 471 

Norris J.R, Marani G.F, Bonnell J.T, 2000, ApJ, 534, 248 

Norris J.R, 2002, ApJ, 579, 386 

Paciesas W.S. et al., 1999, ApJS, 122(2), 465 

Press W.H., Flannery B.P, Teukolsky S.A., Vetterling W.T, 1993, 
Numerical Recipes in C, Cambridge University Press, 2nd edi- 
tion 

Ramirez-Ruiz E. & Lloyd-Ronning N.M., 2002, New A, 7, 197 



